Decoherence can be useful in quantum walks 
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We present a study of the effects of decoherence in the operation of a discrete quantum walk on 
a line, cycle and hypercube. We find high sensitivity to decoherence, increasing with the number 
of steps in the walk, as the particle is becoming more delocalized with each step. However, the 
effect of a small amount of decoherence is to enhance the properties of the quantum walk that are 
desirable for the development of quantum algorithms. Specifically, we observe a highly uniform 
distribution on the line, a very fast mixing time on the cycle, and more reliable hitting times across 
the hypercube. 
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OVERVIEW 



There is currently great interest within the quantum 
information community in quantum versions of random 
walks, because of the possibility they may produce new, 
powerful types of quantum algorithms. All the quantum 
algorithms known until very recently are essentially based 
on application of the quantum fourier transform, includ- 
ing Grover's search [|| and Shor's factoring ||^, which 
provide quadratic and exponential speed up respectively 
over the best known classicl algorithms. Other types of 
quantum computation, such as quantum adiabatic com- 
putation Q, have not yet been shown to provide expo- 
nential speed up over classical methods. Some of the 
most powerful known classical algorithms are based on 
classical random walks, so it is a natural question to ask 
whether there are quantum counterparts that can do even 
better. For example, a random walk on a general graph 
can be used to address hard problems such as approx- 
imating the permanent, fc-SAT and graph connectivity 
i i f. 

Before attempting to create quantum algorithms from 
quantum walks, it is first useful to study their properties 
and dynamics on simpler structures. Several quantum 
analogues of a classical random walk on discrete lattices 
or graphs have been proposed. These include discrete 
time walks both with and without a quantum coin |j , 
and a continuous time walk The relationship be- 
tween the continuous time quantum walk and discrete 
time quantum walks is not fully understood. For the 
cases where they have been studied on the same graph, 
they give essentially the same results, e. g., Quan- 
tum walks on the infinite line, the cycle, and the hyper- 
cube have all been solved analytically |ll|, O, and 
some bounds are known for more general graphs . On 
the line and cycle, for most quantities of interest, such 
as the standard deviation and the mixing time, there is 
a quadratic speed up over the classical walk. Kempe 



p3| recently proved that the hitting time to the oppo- 
site corner of a hypercube shows an exponential speed 
up (a possibility also found numerically by Yamasaki et 
al. ||l4| ). Making an algorithm out of a quantum walk re- 
quires significant further work. Two have been proposed 
very recently, Childs et al. |15 use a continuous time 
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quantum walk to traverse a special graph exponentially 
faster than a classical algorithm, and Shenvi et al. [l^ 
prove that a coined quantum walk can match Grover's 
algorithm searching an unstructured database. 

Quantum walks themselves (both discrete and contin- 
uous time) can be implemented efficiently on a quantum 
computer ISl i. e., it is not necessary to pro- 
vide a physical implementation of a quantum walk to 
base an algorithm on them. Nonetheless, quantum walks 
are interesting in their own right as physical systems in 
which a precise level of coherent control can be demon- 
strated. Several direct implementations of discrete walks 
have been proposed, all with a quantum coin: a walk in 
the vibration modes of a trapped ion |2^, in the phase 
of the field in a cavity containing an atom pT| , and with 
an atom hopping between traps in an optical lattice [ p2[ . 
Other than for these physical implementations, and re- 
cent work on coin decoherence by Brun et al. |2^, p5[ , 
the effect of decoherence in quantum walks has not pre- 
viously been studied in any detail. 

The key observation of this paper is that, in small 
doses, rather than degrading the quantum features, de- 
coherence in a coined quantum walk can enhance the 
desirable quantum speed up, even though overall, quan- 
tum walks (consisting as they do of extremely delocalized 
quantum particles) are highly sensitive to the effects of 
decoherence. This is very encouraging for the prospects 
of using quantum walks as the basis of powerful quantum 
algorithms. 

The paper is organized as follows. First we review one 
of the simplest examples of a quantum walk, the coined 
walk on a discrete line, and describe the properties of 
the perfect quantum walk. Then we present our results 
showing the effects of decoherence in the quantum walk 
on a line, on a cycle, on a hypercube, and on the "glued 
trees" graph of ref . pa . 
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II. COINED QUANTUM WALK 

We consider only coined quantum walks on discrete 
lattices in this paper. Since the classical random walk 
requires a source of randomness (coin-toss) in the dy- 
namics, introducing a quantum coin is a natural way to 
proceed. For the walk on an infinite line, the total Hilbert 
space is H = C'^iS)H°°, where H°° has supposrt on a; G Z. 
We label the coin states by {| —1), | -1-1)} for "move left" 
and "move right" respectively, and those of the particle 
by |a;) for position on the line. We write the tensor prod- 
uct states as \a,x), where a £ {il} is the state of the 
coin. The unitary operator describing a single step of 
the walk is 



U = S-{H(g)l) 



(1) 



where S is the conditional shift operator S\a,x) 
|a, X + a), and H is the Hadamard operator 



{H ig) 1)\a,x) = {a\a,x) + 



(2) 
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FIG. 1; Distribution of the particle position for a quantum 
walk on a line after T = 200 time steps. Pure quantum (dot- 
ted), fully classical (dashed), and decoherence at rate shown 
on part of system indicated by key (solid). Uniform distribu- 
tion between —T/y/2 crosses) also shown. 



acting on the coin states only, it plays the role of a "coin 
toss" . It can be shown for unbiased quantum walks on a 
line in which the "coin toss" prepares an equal superposi- 
tion of states, all such operators are essentially equivalent 
to the Hadamard In higher dimensions, when 

the particle has more than two choices of direction at each 
step, the choice of unbiased coin-toss operator becomes 
correspondingly richer |27, 28|. Unlike the classical case, 
where each coin toss is independent of the previous coin 
tosses, in the quantum case, unitarity of the evolution 
and hence reversibility implies the initial state of the coin 
has observable consequences at all later times. Specifi- 
cally, in the following (unless stated otherwise) we will 
choose the initial state to be (| -|-l)-|-i| —1)) /a/2 which re- 
sults in a symmetric probability distribution ^ . 
The dynamics of the quantum walk thus consist of re- 
peated application of the operator U to the particle and 
coin, resulting in a spreading out on the line, with inter- 
ference causing the quantum speed up. 

The quantum walk on a line has been solved exactly 
| p2[ using both real space (path counting) and Fourier 
space methods. The solutions are complicated, mainly 
due to the "parity" property, i. e., the solutions must 
have support only on even(odd)-numbered lattice sites 
at even(odd) times. The shape of the probability dis- 
tribution for the particle position consists of a nearly 
flat region around the center with the same width as 
the classical binomial distribution, and oscillating peaks 
out towards x = ±T/^/2. Both quantum and clas- 
sical distributions are shown in fig. |l|, calculated nu- 
merically for T = 200. The moments can be calcu- 
lated, for a walk starting at the origin, (|a;|) = T/2 and 
{x^) = (1 - 1/V2)T^ = ct2(T). The standard deviation 
(from the origin) a(T) is thus linear in T, in contrast to 
\/T for the classical walk. 



III. DECOHERENCE IN A QUANTUM WALK 
ON A LINE 



In order to model decoherence in this system, we write 
the coin-particle dynamics in terms of a density matrix 
p{t) that evolves according to 



p{t + l) = (1 - p)U ■ p{t) ■ 

+ pJ2^^■U■p{t)■U^ 



(3) 



Here Pi is a projection that represents the action of the 
decoherence and p is the probability of a decoherence 
event happening per time step. We took equation (|^) 
and evolved it numerically for various choices of . Mo- 
tivated by the likely form of experimental errors, we also 
modelled an imperfect Hadamard by applying a Gaussian 
spread of standard deviation y^7r/4 about the perfect 
value of 7r/2 implicit eq. (|^), compare ||2^. An imperfect 
shift on the particle has been studied in |22|| . In each case 
we find the same general form for the decay of <Jp{T) from 
the quantum to the classical value, with small differences 
in the rates, as shown in fig. ||. The slope of Op{T) is 
finite as p ^ and zero at p = 1. We calculated crpiT) 
analytically (details in ref. [g9|) for pT <C 1 and T ^ 1 
for the case where Pi is the projector onto the prefered 
basis {|a,a;)} (decohering both particle and coin), 



ap{T) a{T) 1 - 



6V2 



0{p) 



(4) 



This compares well with simulation data, once a second 
order correction for a{T) = (1 - 1/^2) ^/^(r - l/T) is 
taken into account. The first order dependence is thus 
proportional to pT, so the sensitivity to decoherence 
grows linearly in T for a given decoherence rate p. 
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FIG. 2: Standard deviation o-p{T) of tiie particle on a line for 
different models of decoherence, for T = 100 time steps. 

The quantum walk on infinite regular lattices of higher 
dimension shows a similar decoherence profile, as might 
be expected, since the standard deviation of a classical 
random walk is Vt, independent of dimension. 

IV. DISTRIBUTION SHAPE ON A LINE 

The transition from quantum to classical in the stan- 
dard deviation of the particle position shown in fig. ^ is 
qualitatively the same for all types of decoherence ex- 
amined. However, there are interesting differences in the 
shape of the distribution of the particle position. The de- 
coherence rate that gives the closest to uniform distribu- 
tion has been selected and plotted in fig. |], along with the 
pure quantum and classical distributions for comparison. 
When the particle position is subject to decoherence that 
tends to localize the particle in the standard basis, this 
produces a highly uniform distribution between ±T/ y/2 
for a particular choice of p. The optimal decoherence rate 
Pu can be obtained by calculating the total variational 
distance between the actual and uniform distributions, 

iy{p,T) = ||P(x,p,r)-P„(T)||tv 

^ Y,\Pi^,P,T)-Pu{T)\, (5) 

where P{x,p,T) is the probability of finding the particle 
at position x after T time steps, regardless of coin state, 
and P„(T) = V2/T for -T/V2 < x < T/V^ and zero 
otherwise. Figure ^ shows i'{p, T) for T — 200 with de- 
coherence applied to the coin, particle and both at once. 
Decoherence on both particle and coin produces the best 
uniform distribution, but at a cost of a lower decoher- 
ence rate, and sharper minimum, i. e., greater sensitivity 
to the value of p. Walks with decoherence only on the 
particle can tolerate more variation in the exact decoher- 
ence rate, while not achieving such a good uniform distri- 
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FIG. 3: The total variational distance of the particle position 
distribution from the uniform distribution for T=200 and de- 
coherence on coin (dotted), particle (dashed) or both (solid). 



bution. The optimum decoherence rate depends on the 
number of steps in the walk, we determined numerically 
that PuT ~ 2.6 for decoherence on both and PuT ~ 5 
for decoherence on the particle only. These differences 
in the quality of the uniform distribution are indepen- 
dent of p and T, and provide an order of magnitude (0.6 
down to 0.06) improvement in v over the pure quantum 
value. Decoherence just on the coin doesn't enhance the 
uniformity of the distribution, as fig. ^ shows, there is 
a cusp at X = 0. However, for finite T, there is still a 
useful window within which even coin only decoherence 
does not significantly degrade the linear spreading of the 
walk. 



V. DECOHERENCE ON A CYCLE 

If, instead of an infinite line, the quantum particle is 
allowed to walk on a cycle of size N points, the appro- 
priate quantity to measure the progress of the walk is 
the mixing time. We immediately have to modify the 
classical definition, because, unlike the classical random 
walk which mixes to a uniform distribution in the the 
long time limit, a unitary process such as that evolving 
the quantum walk does not mix to any limit at large 
times. We can instead define a time averaged particle 
distribution 

P(a;,p,T) = -^P(a;,p,t) (6) 

which does always mix for large enough times T |Tl| |. It 
is easy to sample from the distribution P{x,p,T): run 
the walk for some randomly chosen number of steps < 
t < T and measure the particle position at that time t. 
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The mixing time is then defined as 



{r|vt 



>T: \\P{x,p,t)-Pu\\t. <e 



} (7) 



where P„ is the limiting (uniform) distribution over the 
cycle. The mixing time quantifies how long it takes for 
the time-averaged probability distribution of the particle 
position to reach its limiting value within a margin of 
small parameter e > 0. 

The walk on a cycle is the simplest example of a walk 
on the Cayley graph of an Abelian group, and was his- 
torically the first to be treated analytically [0. The 
dynamics of the walk on the cycle are the same as for the 
walk on a line, with the particle position taken mod(iV). 
Aharonov et al. proved an upper bound for the mix- 
ing time of O(A^logA^). The limiting distribution thus 
obtained depends on the choice of coin operator in 
sharp contrast to the classical walk, which always mixes 
to a uniform distribution. For the Hadamard coin used 
here, the odd-iV cycle mixes to the uniform distribution, 
but the even-A^ cycle does not. 

We numerically evaluated the mixing times for walks 
on cycles of sizes up to A^ ~ 80, both for pure states, 
and in the presence of the same types of decoherence 
as described in the previous section for the walk on a 
line. For odd-A^ cycles with no decoherence, we find that 
Aff ~ N/e as compared to the upper bound of ~ 
N\ogN/e in ref. While we believe that the linear 

scaling with TV is the correct result, obtaining a tighter 
bound analytically is a tough task, because the time aver- 
aged probability distribution P(x,p,T) is a rather fluc- 
tuating quantity, especially for p = 0, as is illustrated 
in fig. ^. In this figure, the quantity \\P{x,p,t) — P„||tv 
from eq. (0) is plotted against time. The time at which 
the curves last cross the horizontal line at e = 0.01 is the 
mixing time as plotted in fig. ^. A different choice of e 
thus causes a jump in the value of if it happens to 
touch the next peak in \ \P(x,p,t) — Pu||tv 



A. Even-A'^ cycles with decoherence 
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FIG. 4: Difference b etween the time averaged probability dis- 
tribution P{x,p, T) and the uniform distribution expressed as 
the total variational distance, eqs. for the case with 

decoherence on both coin and particle. Both axes logarith- 
mic. The value of e used in fig. ^ is shown as a horizontal 
line. Top three solid lines labeled on right are for N = 22 
with p = (A), p = 0.001 (B), p = 0.02 (C). Lower four 
lines are for A = 21 with p = (solid), p = 0.002 (dotted), 
p — 0.01 (short dashed), p = 0.02 (long dashed). 
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Under the action of a small amount of decoherence, the 
mixing time becomes shorter for all cases, typical results 
are shown in fig. ^. Also, decoherence causes the even-A^ 
cycle to mix to the uniform distribution. The asymptotes 
in fig. p^for A^ even and decoherence on the coin only, for 
p < 2/N, are well fitted by epM^ ~ A^/4 for A^ divisible 
by 2, and epM^ ~ N/16 for A^ divisible by 4. For larger 
p, the mixing time tends to the classical value of A^^/16e 
[note this is not log(l/e) because we are calculating the 
average mixing time, eq. (^]. Although for A^ divisible 
by 4, the (coin-decohered) mixing time shows a minimum 
below the classical value at p ~ 2/N, this mixing time 
is > N'^/32e, i. e., still quadratic in A^. Thus, although 
noise on the coin causes the even-A^ cycle to mix to the 
uniform distribution, it does not produce a significant 
speed up over the classical random walk. 



FIG. 5: Numerical data for mixing times on cycles of size 
iV = 29 and A = 30 (□), for coin (dotted), particle (dashed) 
and both (solid) subject to decoherence, using e — 0.01. Also 
A = 28 ( X ) for coin. Both axes logarithmic. 



For decoherence on the particle position, with p < 
16/N^, epM^ ~ 1/(^/2 - 1) for A^ divisible by 2, and 
epM, ~ l/{N/4 + 3) for A^ divisible by 4. At p ~ 16 /N^, 
there is a minimum in the mixing time at a value roughly 
equal to the (A^ ± l)-cycle pure quantum mixing time, 

Afi™'"'' ~ aN/e (with a a constant of order unity). The 
top three lines in fig. ^ show how decoherence pulls an 
even-N cycle down to mix to the uniform distribution at 
the same rate as the neighboring odd-A^ cycle. Decoher- 
ence on the particle position thus causes the even-A^ cycle 
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to mix to uniform in linear time for a suitable choice of 
decoherence rate ^ independent of e. 



B. Odd- AT cycles with decoherence 

For all types of decoherence, the odd-A^ cycle shows 
a minimum mixing time at a position somewhat earlier 
than the cvcn-iV cycle, roughly p = but because 

of the oscillatory nature of P{x,p,T), the exact behav- 
ior is not a smooth function of p or e. As decoherence 
on the particle (or both) increases, the oscillations in 
P{x,p, T) are damped out. The lower set of lines on fig. ^ 
shows how both these features affect the mixing time. At 
p ~ the mixing time passes smoothly through an 

inflexion and from then on behaves in a quantitatively 
similar manner to the adjacent-sized even-A'^ cycles, in- 
cluding scaling as M^™™^ ~ aN/e at the inflexion. Thus 
for < p < 16/A^^ there is a region where the mixing 
time stays linear in A^. Our overall conclusion is thus the 
same as for the walk on a line, there is a useful window 
within which decoherence enhances rather than degrades 
the quantum features of the walk. 



VI. HYPERCUBE DECOHERENCE 

The hypercube (boolean A^-cube) is also the Cayley 
graph of a group, and provides a step on the way to quan- 
tum walks on more complex structures. Since there are N 
edges joining at each vertex, we need an Af-dimensional 
coin to choose between the possible paths at each step. 
This opens up a correspondingly larger range of possible 
unitary operations to use for the coin toss, but a sensible 
choice is one which respects the symmetry of the under- 
lying graph, and is as far from the identity as possible. 
For the hypercube, this is the Grover operator, whose 
elements expressed as an N x N matrix are defined as 
Gjk = 2/A^ — djk- We also choose a symmetric starting 
state for the coin with equal weights for all possible di- 
rections from the chosen starting node. The high degree 
of symmetry in the hypercube allows the quantum walk 
to be mapped to a walk on a line with a different coin 
toss operation at each point . While the mixing time 
for a quantum walk on a hypercube is somewhat worse 
than a classical random walk, Kempe proved that 
the hitting time to the opposite corner is polynomial, an 
exponential speed up over the classical walk. 

Kempe discusses two types of hitting times, one-shot, 
where a measurement is made after a pre-determined 
number of steps, and concurrent, where the desired loca- 
tion is monitored continuously to see if the particle has 
arrived. In each case, the key parameter is the probability 
P/i of finding the particle at the chosen location. We cal- 
culated Ph numerically and found that all forms of deco- 
herence have a similar effect on Ph, see fig. ^ reducing the 
peaks and smoothing out the troughs. For the one-shot 




time step 



FIG. 6: Hitting probability on a 9-dimensional hypercube for 
one-shot (left) and concurrent (right), perfect walk (circles), 
with p = 0.05 (dotted), p = 0.1 ~ 1/9 (solid). Classical 
hitting probability barely visible (dashed). 



hitting time this is useful, raising Ph in the trough to well 
above the classical value, so it is no longer necessary to 
know exactly when to measure. For p < 1/N, the height 
of the first peak scales as Ph{p) = Ph{0) exp{—{N + a)p}, 
where < a < 2 depending on whether coin, particle 
or both are subject to decoherence. An exponential de- 
crease in the presence of decoherence sounds about as bad 
as it could reasonably be, and for long times, of course, 
decoherence reduces the walk to classical behavior. How- 
ever, the hitting times are short, only ~ N-k /2 steps, and 
p ~1/N only lowers Ph by a factor of I/e. For algorith- 
mic purposes this is insignificant, only a factor of order 
unity and thus still exponentially better than classical. 
(Standard amplification techniques can be used to bring 
the hitting probability as close to one as desired.) Note 
also, that the size of the graph (measured in number of 
nodes) is exponential in N , so the decoherence has only 
a linear effect measured in terms of the size of the graph. 

Continuous monitoring of the target location as in the 
concurrent hitting time is already a sort of controlled 
decoherence, note that the height of the initial peak at 
~ Ntt/2 steps is only about 2/3 that of the one-shot 
hitting probability. No extra features are produced by 
the addition of unselective decoherence, but there is still 
a range of < p < \/N within which the quantum 
speed up is preserved. Note that in both the one-shot 
and concurrent cases, p ~ l/N \s & critical damping 
rate, smoothing out the second peak (shown at around 
40 (ee 3A/7r/2) steps in fig. |). 

VII. "GLUED TREES" WALK 

The graph used in Childs et al. is similar to the 
hypercube in that it can also be mapped to a walk on 
a line. It consists of two identical binary trees of depth 
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A'', glued at their branches by a randomly distributed 
set of edges (two per node so the degree of the graph 
is three except for the roots of the trees). They use a 
continuous time walk for their algorithm, but a discrete 
quantum walk using a three dimensional Grover coin op- 
erator has essentially similar properties |2^. The 
effects of decoherence are similar to the hypercube. The 
peak in the probability of finding the particle at the exit 
node occurs after 27V + 1 steps, and is around 0.6 for 
a pure walk. Decoherence reduces this peak exponen- 
tially (in N) while spreading out the range of time steps 
over which the probability is significantly (exponentially) 
larger than the classical value. 

VIII. CONCLUSIONS 

One of the generic ways in which classical random 
walks are applied to algorithms is to sample an exponen- 
tially large problem space to estimate statistical proper- 
ties of the system. Fast mixing times and uniform sam- 
pling are necessary properties of the random walk process 
for it to perform efficiently. This is the sense in which 
we propose that a small amount of decoherence in a walk 
on the line or the cycle is beneficial, producing more uni- 
form distributions (line) and faster mixing to a uniform 
distribution (cycle). 

On the hypercube and "glued trees" walks, the key 
quantum property exploited by the examples in the liter- 
ature is the opposite of a uniform distribution, the ability 
of a quantum walk to continue its forward march through 
the graph despite the many possible "wrong turns" a clas- 
sical random walk gets lost in. For the examples given, 
it is easy to determine a priori the best time to check 



for the particle having found the target node. But for 
more complex graphs this may not be easy to calculate, 
so having a wider window of opportunity to successfully 
detect the particle could be an advantage. Note that the 
concurrent hitting time is doing exactly this in a more 
precise way by monitoring the target node at every time 
step. For these type of walks, we don't claim that de- 
coherence gives a major advantage, only that it is not 
detrimental for small decoherence rates. 

We have also found numerically that the walk on a 
cycle mixes in linear time (compared to the upper bound 
of 0{N log A^) proved in ||ll|), and show why it is hard to 
prove a tighter bound analytically because the quantities 
involved are highly fluctuating. 

Overall, these results are quite promising for the de- 
velopment of further quantum algorithms, and for the 
practical implementation of quantum computing. They 
are also encouraging for the experimental implementa- 
tion of quantum walks as proposed in rcfs. |po| , p2[ . 
For the modest number of steps in the proposals, our re- 
sults suggest a very reasonable range in which quantum 
effects can be observed experimentally. Finally, it is also 
an intriguingly counter- intuitive result in its own right, 
worthy of further study in the context of the transition 
from quantum to classical. 
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